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Abstract 

The nonlinear effect of gluon depletion in the collision of large nuclei can be large. It 
is due to multiple scatterings among comoving partons initiated by primary scattering 
of partons in the colliding nuclei. The effect can give rise to substantial suppression of 
J/U production in very large nuclei, even if the linear depletion effect is insignificant 
for the collisions of nuclei of smaller sizes. This mechanism offers a natural explanation 
of the enhanced suppression in the Pb-Pb data recently observed by NASO. 


In a previous paper [Q] we have examined the issues involved in ascribing some aspect 
of the phenomenon oi J/'ip suppression in heavy-ion collisions P] to the depletion of gluons 
prior to the hard subprocess of cc production. What we have found is that the data on the 
survival probability S without the points from Pb-Pb collisions [§, by themselves, cannot 
distinguish whether the suppression is due to gluon depletion or hadronic-nuclear absorption. 
That is, both mechanisms contribute to an exponential dependence of S on the effective path 
length L (or on \ogAB). We now consider the enhanced suppression in the Pb-Pb data of 
NA50 [0 and show how gluon depletion can naturally account for it. Furthermore, it is 
possible for that to happen even if the “normal” suppression in the lighter-ion data is due 
mainly to the absorption mechanism with negligible depletion effect. 

Many suggestions have been advanced to account for the enhancement oi J/'ip suppression 
observed in the Pb-Pb collision data i-in. They all refer to the absorption processes after 
the production of the cc state. Our suggestion is concerned with the depletion of gluons 
before the gg —>• cc subprocess. The basic idea is rather intuitive and can be described 
qualitatively before we go into the details. Consider a row of nucleons in nucleus A colliding 
with another row in nucleus B, and suppose that the uaAi one from the front of the former 
(call it a) collides with the n^th one in the latter (call it h) in a hard process creating cc. 
The gluon depletion mechanism discussed in |Q takes into account of the loss of gluons in 
a (due primarily io g ^ qq) as it goes through B until gg cc occurs with a gluon in 6; 
similarly, the gluons in b are depleted as b traverses A. We shall refer to this process as linear 
depletion for reasons that will become clear below. What we now want to emphasize is that 
a nonlinear depletion process may be even more important. Such a process is due to the 
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interaction of the gluons in a with the slower partons liberated from the ha — 1 forerunners 
in A broken by earlier interactions, and likewise b with the partons of the — 1 forerunners 
in B. In an imperfect, yet helpful, analogy one may think of a multicar accident on a busy, 
foggy highway and recognize that most of the collisions are between cars originally going in 
the same direction. 

Let us hrst summarize the essence of the linear effect considered in . Our notation will 
follow that of |l|, but abbreviated. The probability that a nucleon in A makes ui collisions 
in B before the hard subprocess is 


TTvi = —^ ns = CTinTg 

z/i! 


( 1 ) 


where is the path length that is traversed in B before the cc production and is dependent 
on the impact parameter bs and longitudinal position zb, both being suppressed. If the 
depletion factor per collision at hxed momentum fractions xi and X 2 is D {D = 1 for no 
depletion), then the suppression factor at hxed b and z in A and B is 

Tab = Z! = exp [-(1 - D){nA + ub)] , (2) 

where is dehned as in (|^), but with B replaced by A. It is the simple sum, + n^, in 
(0) that leads us to call the effect linear. The exponential behavior of T^^^ is what generates, 
after integration over b and z, the approximate exponential dependence of S' on L that 
is indistinguishable from the Gerschel-Hufner formula [^, derived from purely absorptive 
consideration. 

The nonlinear effect that we now describe arises from the interactions with the forerun¬ 
ners. At the partonic level the linear effect is due to the primary interactions of gluons 
in nucleons going in opposite directions, while the nonlinear effect is due to the secondary, 
tertiary, etc. interactions among partons moving in the same direction, initiated by primary 
interactions. The rapidity separation. Ay, between the participants of the primary inter¬ 
action is large because they belong to the nuclei A and B separately. On the other hand. 
Ay between the partons involved in the secondary (or tertiary, etc.) interactions is small 
because they belong to the same nucleus. Ordinarily, in an unperturbed nucleus or in deep 
inelastic scattering of a nucleus, those partons in different nucleons do not interact except in 
the context of nuclear binding and shadowing. However, if a primary interaction has taken 
place between two colliding nuclei, the scattered parton in A, whether at large or small angle, 
can interact with a parton coming from behind in the same or neighboring rows. Since they 
are comovers, their interaction can be much stronger than the primary interaction, a prop¬ 
erty that is consistent with the general notion of strong interaction in soft processes being 
short-ranged (in rapidity). Thus even if the linear depletion effect is small, the nonlinear 
effect need not be. 

If we consider an xn^ matrix, representing the possible pairings of ua and ub nucleons 
in collisions, the last row and last column contribute to the linear depletion effect. [Their 
sum Ua + tib — 1 appears as ua + ns in (@) in compensation for the fact that the hrst 
collision of a nucleon with a row of nucleons is the normal pp collision, whose cross section 
is larger than those of the subsequent collisions that involve the broken nucleon propagating 
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downstream. To elaborate on this point is too mnch of a digression that is not germane to 
the following discussion.] The remaining part of the matrix having {ua — 1)(^b — 1) pairings 
contributes to the quadratic depletion effect due to multiple parton scatterings. Let us dehne 

u'a = {ua - 1) Q{nA - 1) , (3) 


and similarly for n'^. Then, assuming A < B, the average number of collisions that the 
forerunners of a in A make with the forerunners of 6 in i?, producing comoving partons that 
can interact with the partons of a, is n'^n'^ — n'j^/2] that for producing comoving partons 
with the ones in b is j2. This way of partitioning the n’An’^ pairings can be visualized in 
the forward light-cone of AB collision, where the former lie on the A side of the interaction 
region, while the latter lie on the B side. The precise method of partitioning is unimportant, 
as will become evident presently. 

The probabilities that a and b can interact and times with their respective forerun¬ 
ners are 


— 


U[\ 
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If D' is the effective gluon depletion factor for each of those interactions, then the corre¬ 
sponding suppression factor, analogous to (||), is 


^ AB — 


E TT I 7T I 
^1 ^2 


D' 


exp [-(1 - D') n'^n's] 


( 6 ) 


We refer to this as the quadratic depletion effect, since it is that appears in the 

exponent, as opposed to ua + ub in (H). As it is in (H), the dependences on 6^, za, and 
Zb have been suppressed in (H). 

The combined suppression factor due to both linear and quadratic depletion as well as 
absorption is now 

P = eyi]i [-{I-D){nA + nB) - {I-D')n'xn'B-aairX+ T^)] , (7) 

where cXa is the absorption cross section and is the path length in A traversed by the 
J/^l] system. Exhibiting the b and 2 ; dependences, we have [|[] 

T^ = {l-\) Po{La±za) , LA = {R\-s^f/^ , ( 8 ) 

and similarly for with bA = s'and bB = b — s. The average overall suppression factor 
(more precisely, survival probability) is 

= Nxh I d% J (fs dzA f_" dzB P , (9) 
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where Nab is the same integral as in (||) but with P replaced by 1. 

To see how depends on A and B, let us examine the parameters in the formula. 
Without the quadratic depletion terms in (|^), we have 

P, = P{D' = 1) = exp l-aa (T^ + Tg) - (Ti + T+)] , (10) 

where ad = cTin(l — D). As pointed out in |I|, (pT]) exhibits the symmetry between the 
depletion effect before the formation of J/ijj and the absorption effect afterwards. That is 
why the exponential dependence of the empirical on the effective length L (or \ogAB) 
cannot distinguish the two effects. So long as the combined cross section ac = aa + ad is 
around 7 mb, the heavy-ion data, excluding the Pb-Pb collisions, can be htted by any ratio 
r] = ad/aa- Now, we consider the contribution from the quadratic depletion term in (|^) only, 
giving 


P 2 = P{D = 1, (ja = 0) = exp [- T {riA - 1) {ub - 1) <d{nA - 1) 6(^5 - 1)] 


( 11 ) 


where r = 1 — D', a parametrization, like ad, having the more proper sense of depletion in 
that r = 0 means no depletion. There are a number of features of (|Tl|) worth noting. 

(a) While the discussions in the introduction and in the paragraph containing (^) regard 
Ha and hb as integers for the sake of ease in describing the nonlinear depletion mechanism, 
they can in reality have any positive value by virtue of their dehnitions, nA,B = crinTl^. 
That is why the step functions in and ([TI|) are important to ensure that the participants 
of the process, ti'a and n'^, are nonnegative. As a consequence there is a threshold effect, 
i.e., A and B must be large enough for the mechanism to be operative. 

(b) The inelastic cross section cxin is relevant in the determination of the position of the 
threshold. It is not the afl for pp collision because, except for the hrst collisions on the 
front sides of the nuclei, most of the collisions are between broken nucleons 
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consist mainly of the parton fluxes that propagate downstream after the bound nucleons 
are broken by the hrst collisions. Uin is an effective cross section for the collision of such 
broken nucleons, and there exist no reliable estimates for its value. Using p' to denote broken 
nucleon, and taking cr-Jf ~ 30 mb, it is not unreasonable to consider afj’ ~ 20-25 mb, and 
~ ~ 13-21 mb. We shall adopt Uin ~ 15-25 mb as typical values. 

(c) The quadratic depletion parameter r can be substantially different from zero, even if 
the linear effect measured by ad is zero, since, as discussed earlier, the interaction between 
partons with small rapidity separation can be much greater than that between partons with 
large Ny. Since the determination of r from hrst principles is difficult, we shall use it 
as a free parameter in the following. It should be noted that even if D' = 0, i.e., total 
depletion per collision, r attains its maximum value 1 , so ([Tl|) does not give P 2 = 0 . That 
is because the Poissonian huctuations in (^) and (|^) allow for z/( = 1/3 = 0, which result in a 
nonvanishing probability for the passage of the gluon huxes with minimal inhuence by the 
depletion mechanism. 

To gain some further insight in the quadratic depletion ehect, let us compute taking 
only P 2 into account, i.e., by substituting (0) alone into (^). Using the integration procedure 
developed in |^, we obtain the results shown in Fig. 1, where Uin is set at 15, 20 and 25 
mb; the shaded regions are bounded by r = 0.5 from above (for illustrative purpose) and 
r = 1.0 from below. Evidently, the threshold for the quadratic depletion effect is higher 
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at smaller a-m, since there would be less participants for the multiscattering subprocesses 
unless A is higher. Furthermore, even at maximum depletion (r = 1) there is still a residual 
rate of J/-^ production because of the aforementioned probability of gluon passage without 
depletion. We note that, although the parameters Uin and r are not empirically familiar, 
their values used in Fig. 1 are sensible estimates, so the suppression effect revealed is a 
natural consequence of a physical process that is not contrived to explain the data. 

For a comparison with the data [P, Q we include both Pi and P 2 in (P) and calculate 
the overall suppression factor. Since, as found in |^, the combined effect of absorption and 
linear depletion is insensitive to the ratio rj = ad/ua, we choose the uncontroversial values: 

+ o'rf = 7 mb, and rj 0.1. For quadratic depletion effect we use (Tin = 20 mb and 
r = 0.5-1.0. The result is shown by the triangles in Fig. 2. The agreement with the data 
[4] is evidently very good. Of course, if there exists enhanced nuclear, hadronic or plasma 
absorption at high AB, it can be accommodated by reducing the value of r. What is shown 
here is that the quadratic gluon depletion effect by itself is able to account for the enhanced 
suppression in the Pb-Pb data. 

A concomitant phenomenon associated with quadratic gluon depletion is the suppres¬ 
sion of back-to-back DD production in Pb-Pb collision, but not in AB collisions where A is 
smaller. Photon production would not necessarily be suppressed, since the quarks produced 
by gluon conversion can carry on the 7 -producing subprocess without inhibition. Dilepton 
production may or may not be enhanced, depending on whether the extra quarks and anti¬ 
quarks are produced inside or outside the interaction region. It is therefore important that 
all those signatures should be examined experimentally in the collisions of very heavy ions. 

Whether or not the nonlinear gluon depletion process can wholly or partially account 
for the enhanced J/'ip suppression phenomenon, what we have discovered here is that there 
is a whole class of parton interactions whose role in heavy-ion collisions has hitherto been 
overlooked, but they are of crucial importance to any process whose rate depends on the 
magnitude of the gluon flux available in large nuclei. 

This work was supported, in part, by the U.S.-Slovakia Science and Technology Program, 
National Science Foundation under Grant No. INT-9319091 and by the U. S. Department 
of Energy under Grant No. DE-FG03-96ER40972. 
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Figure Captions 

Fig. 1 The snppression factor abbreviated as S', is plotted against AB, when only 
the qnadratic depletion effect is taken into acconnt. The shaded regions are for r having 
valnes between 0.5 (npper bonndaries) and 1.0 (lower bonndaries). Three valnes of crj„ 
are nsed, as indicated. 

Fig. 2 The snppression factor S'jy^, abbreviated as S', is plotted against AR, when both 
the nsnal linear (mainly absorption) effect and the qnadratic depletion effect are taken 
into acconnt. The data are from [4]. Typical valnes for the parameters in the theoretical 
calcnlations have been nsed. 
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